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ABSTRACT 

In this Paper we report on further observations of the third and fourth kilohertz 
quasi-periodic oscillations (QPOs) in the power spectrum of the low-mass X-ray bi- 
nary (LMXB) 4U 1636-53. These kilohertz QPOs arc sidebands to the lower kilohertz 
QPO. The upper sideband has a frequency 55.5±1.7 Hz larger than that of the contem- 
poraneously measured lower kilohertz QPO. Such a sideband has now been measured 
at a significance > 6<r in the power spectra of three neutron star LMXBs (4U 1636-53, 
4U 1728-34, and 4U 1608-52). We also confirm the presence of a sideband at a fre- 
quency ^55 Hz less than the frequency of the lower kilohertz QPO. The lower sideband 
is detected at a 3.5cr level, only when the lower kilohertz QPO frequency is between 
800 and 850 Hz. In that frequency interval the sidebands are consistent with being 
symmetric around the lower kHz QPO frequency. The upper limit to the rms ampli- 
tude of the lower sideband is significant lower than that of the upper sideband for lower 
kilohertz QPO frequencies >850 Hz. Symmetric sidebands are unique to 4U 1636-53. 
This might be explained by the fact that lower kilohertz QPO frequencies as high as 
800-850 Hz are rare for 4U 1728-34 and 4U 1608-52. Finally, we also measured a low 
frequency QPO at a frequency of ~ 43 Hz when the lower kilohertz QPO frequency 
is between 700-850 Hz. A similar low-frequency QPO is present in the power spectra 
of the other two systems for which a sideband has been observed. We briefly discuss 
the possibility that the sideband is caused by Lense-Thirring precession. 
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1 INTRODUCTION 

One of the important discoveries made with the Rossi X- 
ray Timing Explorer (RXTE) satellite is that of the kilo- 
hertz quasi-periodic oscillations (QPOs) in the power spec- 
tra of ^ 20 accreting neutr on star low-mass X-ray bina- 
ries (see Ivan der Klisllioool for a review). Kilohertz QPOs 
are caused by the motion of matter in regions of spacetime 
within a few kilometres of the surface of accreting neutron 
stars, where strong-field gravity is required to describe such 
motion. They potentially allow one to detect strong-field 
effects and to constrain the neutron star mass-radius rela- 
tion. Although there is as yet no agreement about the precise 
physical mechanism underlying kilohertz QPOs, most mod- 
els agree that the frequency of one of the observed kilohertz 
QPOs reflects the frequency of orbital motion at the inner 
edge of the accretion disc. 



The basic phenomenology of kilohertz QPOs consists of 
two kilohertz QPO peaks, which are separated by Ai^=200- 
360 Hz and which move in frequency by up to 700 Hz in gen- 
era l correlation with mass accretion rate indicators (again 
see lvan der Klisl2000l for a review). The highest observed fre- 
quency is ~ 1330 Hz llvan Straaten et afll200Cft . correspond- 
ing in the case of a 1.4 Mg neutron star to an orbital radius 
as tight as 15 km. Some of these sources also show X-ray 
bursts, in which burst oscillations are seen which last for a 
few seconds and have frequencies between 270 and 620 Hz 
(see IStrohmaver fc Bildstenl 120041 for a review) . These os- 
cillations usually drift by 1-2 Hz but are near the neutro n 
star spin frequency in each source llChakrabartv et"al .1200311 . 
In the millisecond pulsar SAX J1808. 4-3658 Av has been 
found to be equal to half the spin frequ ency, which dis- 
prove s spin-orbit beat frequency models iWiinands et alJ 
2003) and suggests instead that the neutron star spin in- 
duces resonances in the disk flow, perhaps involving the 
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general relativistic epicyclic frequencies lAbra mpwicz et alJ 
l2003l:IWiinands et alJl2003t lLamb fc Millerll2003li ~ 

After extensive searches for additional kilohertz QPO 
peaks at theoretically predic ted frequencies which remaine d 
unsuccesful ( e .g. in Sco X-l: lMendez fc van der KlidEicI) . 
Ijonker et alJ J200Q) discovered sidebands at a frequency 
50-65 Hz above the lower kilohertz QPO in 4U 1608-52, 
4U 1728-34, and 4U 1636-53. Magnetospheric modulation 
of the lower kilohertz QPO, a beat phenomenon taking place 
inside the marginally stable orbit as well as Lense-Thirring 
precession are pos si ble exp l anatio ns for these sidebands 
j Jonker et ai1l2000l) . iPsaltisl i2000T) demonstrated that the 
hydrodynamic disk mode model IPsaltis fc Normanlll999l) 
naturally produces a sideband at a frequency near that ob- 
served plus several other, as yet unobserved peaks; no other 
models have so far been able to explicitly accommodate the 
sideband. The 50-65-Hz frequency difference between the 
lower kilohertz QPO and the sideband frequency, the 'side- 
band separation' Avsb, is different in each source and not 
identical to any of the other QPO frequencies simultaneously 
obser ved in the 10-100 Hz frequency range jjonker et alJ 
|2000D . 

In this Paper we find evidence for symmetric sidebands 
to the lower kilohertz QPO in the atoll source 4U 1636-53. 



2 OBSERVATIONS, ANALYSIS, AND 
RESULTS 

We used all the proportional counter array (PCA) data from 
RXTE observations of 4U 1636-53 available to us at the 
beginning of 2004. Hence, we included the data used by 
Ijonker et alJ i2000l) . The analysis we per formed was the 
same as that done by Ijonker et al.l <|2000F) . Below, a con- 
densed descripti on of the analysis ste ps is given; more details 
can be found in Ijonker et alJ (l2000l) . 

Using 128 s-long segments of high time resolution PCA 
data (122/is resolution), we calculated power spectra up to a 
Nyquist frequency of 4096 Hz in an energy band of 2-60 keV. 
The power spectra were searched for lower kilohertz QPOs 
which are narrow (full-width at half maximum, FWHM, less 
than ~10 Hz) and detected in 128 s at a significance larger 
than 2a. This resulted in a selection of 244 ksec of data. 

The lower kilohertz QPO was tr aced using a dynam ical 
power spectrum (e.g. see plate 1 in iBereer et aflil996ft to 
visualise the time evolution of the QPO frequency. For each 
power spectrum the lower kilohertz QPO peak was fitted in a 
range of 200 Hz centred on the traced frequency using a func- 
tion consisting of a constant plus a Lorentzian. This provides 
us with a lower kilohertz QPO frequency measurement for 
each 128 s power sp ectrum. We us e d the shift-and-add tech- 
nique described bv lMendez et all lll99St) to shift each lower 
kilohertz QPO to a reference frequency. Next, the shifted, 
aligned, power spectra were averaged. The average power 
spectrum was finally fitted in the range 512-2048 Hz so as 
to exclude the edges which are distorted due to the shifting. 
The fit function consisted of a constant to fit the noise and 
three or four Lorentzians representing the QPOs. The fit- 
ted average power spectrum is displayed in Figure Errors 
on the fit parameters were calculated using A^ 2 = 1.0 (la 
single parameter). 

We clearly detected the two kilohertz QPOs and the 
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Figure 1. The average power density spectrum of 4U 1636-53 
showing the lower kilohertz QPO (the peak is off— scale), the up- 
per kilohertz QPO, and the upper sideband to the lower kilohertz 
QPO. Due to the shift-and-add method applied in the analysis 
only the frequency differences are meaningfull (see text and refer- 
ences therein). The drawn line represents the best fit of a function 
consisting of a constant and three Lorentzians. The power in the 
Y— axis is given in Leahy units. 

sideband to the lower kilohertz QPO. The sideband was 
detected at a significance of 6a, confirming our previous 
3a d etection of the sideband in 4U 1636-53 ijonker et alJ 
2000). The frequency difference between the frequency of 
the lower kilohertz QPO and that of the sideband, A^sb = 
55.5 ± 1.7 Hz. The average kilohertz QPO frequency separa- 
tion, (Av = Supper- flower), is 283±5 Hz. The FWHM of the 
sideband, the lower kilohertz QPO and the upper kilohertz 
QPO is 12.7±2.7 Hz, 4.86±0.04 Hz, and 130±11 Hz, re- 
spectively. Due to the fact that we have combined data with 
either 3, 4, or 5 active proportional counter units (PCUs) we 
cannot determine the fractional rms amplitude of the QPO 
peaks from the combined data. Therefore, we selected data 
where only 3, 4, or 5 PCUs were active, we fitted those data 
sets separately. In those separate fits we fixed the FWHM 
and the frequency of the QPOs to the values found using the 
complete data set. The fractional rms amplitude of the side- 
ban d is consistent with what was found before (see Table Q 
and Ijonker et afll2000h . 

In order to search for changes in the sideband sepa- 
ration frequency A^sb we averaged shifted power spectra 
based on the unshifted frequency of the lower kilohertz QPO 
in three frequency bands; a lower kilohertz QPO frequency 
between 700-800 Hz, 800-850 Hz, and above 850 Hz. Inter- 
estingly, the average power spectrum in the lower kilohertz 
QPO frequency interval 800-850 Hz showed evidence for two 
sidebands located symmetrically around the lower kilohertz 
QPO peak (see Figure The results of fits to these aver- 
aged power spectra are given in Tabled 

Addin g more data with respect to the analysis per- 
formed bv Ijonker et alJ l)2000l) . we now find a QPO at low 
frequencies (~43 Hz, see Figure[3J. In order to compare the 
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Table 1. Properties of the lower and upper kilohertz QPO (y\ and U2, respectively), the sidebands to the lower kilohertz QPO, and 
the low frequency QPO. The properties displayed below the last horizontal line have been determined without using the shift-and-add 
technique. The FWHM of the lower sideband has in all cases been fixed to that of the upper sideband. 



All frequencies 700-800 Hz 800-850 Hz 850-920 Hz 



Amount of data (ksec) 


244 


Mean i/i±stan. dev. (Hz) 


830±50 


Ai/ (Hz) 


283±5 


A^SB, lower (Hz) 




A^SB, upper (Hz) 


55.5±1.7 


rms amplitude v\ (%) 


7.54±0.01 


FWHM v x 


4.86±0.04 


rms amplitude V2(%) 


4.03±0.11 


FWHM v 2 


130±11 


rms amplitude u^b, lower (%) 


<0.9 a 


rms amplitude i>sb, upper (%) 


1.3±0.1 


FWHM ^SB.upper 


12.7±2.7 



64 


80 


99 


763±25 


828±15 


877±17 


325±6 


286±9 


249±4 




55.1±2.5 




47.7±3.5 


57.0±2.1 


56±2 


7.94±0.04 


7.86±0.02 


7.12±0.03 


5.4±0.1 


4.45±0.06 


4.90±0.06 


4.5±0.3 


4.0±0.2 


3.4±0.2 


104±15 


120±18 


87±13 


<1.5 a 


1.2±0.2 


<0.8 a 


1.7±0.2 


1.4±0.2 


1.2±0.1 


14±5 


12±4 


10±4 



rms amplitude v lo „ (%) 2.7±0.2 2.3±0.1 <1.< 

FWHM y low (Hz) 17±3 20±4 

Frequency v low (Hz) 42.1±0.8 43.5±1.6 

a 95% upper limits (determined using A\ 2 = 2.71). 
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Figure 2. The average power density spectrum of 4U 1636—53 se- 
lecting only lower kilohertz QPO frequencies between 800-850 Hz. 
The two kilohertz QPOs and the sideband on the high frequency 
side of the lower kilohertz QPO are visible. In addition, a side- 
band on the low frequency side of the lower kilohertz QPO can 
be discerned (3.5cr single trial). Due to the shift-and-add method 
applied in the analysis only the frequency differences are usefull 
(see text and references therein). The drawn line represents the 
best fit of a function consisting of a constant and four Lorentzians. 



properties of this low-frequency QPO with that of the kilo- 
hertz QPOs and the kilohertz QPO separation frequencies 
we fitted the average unshifted power spectrum in the fre- 
quency range from 1/16-256 Hz in the three lower kilohertz 
QPO frequency intervals as shown in Table Q In the fre- 
quency range 850 Hz and above the low frequency QPO was 
undetectable (see Table 0. 
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Figure 3. The low frequency part of the average power den- 
sity spectrum of 4U 1636-53 selecting only lower kilohertz QPO 
frequencies between 700-850 Hz. A low frequency QPO at a fre- 
quency of ~ 43 Hz is visible. The drawn line represents the best 
fit of a function consisting of a power law to fit the underlying red 
noise continuum and a Lorentzian to fit the low frequency QPO. 
We did not apply any shift to produce this power spectrum. 



3 DISCUSSION 

Following the discovery by I.Tonker et alJ {2000) of a new, 
third, kilohertz QPO in the power spectra of the three atoll- 
type low-mass X-ray binaries 4U 1608-52, 4U 1728-34, and 
4U 1636-53 we obtained additional data with the RXTE 
satellite of 4U 1636-53 with the goal to study in detail the 
properties of this third kilohertz QPO, also known as the 
sideband to the lower kilohertz QPO. We selected 4U 1636- 
53 since in that source there is evidence for the presence 
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of two sidebands located symmetrically around the lower 
kilohertz QPO peak frequency. 

We confirm the presence of a sideband at a frequency 
55.5±1.7 Hz higher than that of the lower kilohertz QPO 
frequency in 4U 1636-53 at a 6 a significance level. Further- 
more, the addition of the extra data allowed us to investigate 
the sideband (s) as a function of the lower kilohertz QPO fre- 
quency. We found that when the lower kilohertz QPO has a 
frequency in the range 800-850 Hz there is significant evi- 
dence (3.5<t, single trial) for the presence of a sideband at the 
lower frequency side of the lower kilohertz QPO. Selecting 
lower kilohertz QPO frequencies >850 Hz and when com- 
bining all the data, the rms amplitude of the lower sideband 
is significantly less than that of the upper sideband. For the 
frequency range 700-800 Hz the upper limit is consistent 
with the rms amplitude of the upper sideband. 

In t he framework of the son ic-point and spin-resonance 
model of lLamb fc Mille'rl J2003T) . the sideband is unlikely to 
be caused by Lense-Thirring precession at the radius where 
the spin-resonance occurs (near the orbital radius with a 
Kepler frequency of v sp in/2), since for neutron star spin fre- 
quencies such as those considered here, the Lense-Thirring 
precession frequency at that radius is ~2 Hz, this is much too 
low. It is possible though that the radiation pattern responsi- 
ble for the upper kilohertz QPO is modulated by the Lense- 
Thirring precession at the sonic-point radius, this modula- 
tion would then in turn modulate the formation of the lower 
kilohertz QPO. In such a scenario the upper kilohertz QPO 
should also display sidebands. The fact that these have not 
yet been detected can be explained by the large FWHM of 
the upper kilohertz QPO, which precludes the detection of 
such sidebands. The skewed shape of the upper kilohertz 
QPO in the average power spectrum (see e.g. Fig can 
be explained by the fact that the upper kilohertz QPO rms 
amplitude and Av change as a fu nction of kilohertz QPO 
frequency (e.s. iMendez et alJIlflfld Table Q). 

If Lense-Thirring precession is responsible for the pro- 
duction of the sideband, and if the sideband separation fre- 
quency Avsb, reflects the Lense-Thirring precession fre- 
quency, Avsb would be expected to change as the Kep- 
lerian frequency squared (Avsb oc v^). In Figure |1] we 
plotted Avsb as a function of the lower {left panel) and 
upper kilohertz QPO (right panel), respectively. The side- 
band separation frequency Avsb is consistent with being 
constant when the lower/upper kilohertz QPO frequency 
changed by ~ 100/~ 38 Hz, respectively. The best fitting 
quadratic relation that would be appropriat assuming the 
lower/upper kilohertz QPO reflects the Keplerian motion 
at a preferred radius in the disc and Avsb the Lense- 
Thirring precession frequency at that same radius is also 
shown. The normalisation of the quadratic curve is given 

by z^lt = 8 ^3 (^t i s the Lense-Thirring precession 

frequency, I the moment of inertia of the neutron star, v s 
the spin frequency of the star which is ^ 581 Hz in case 
of 4U 1636-53 from IStrohmaver et al.lll997l. c the speed of 
light, and Mns the mass of the neutron star). The curve 
in the left panel requires 145/m to be 3.05±0.07, that in 
the right panel 0.17±0.01, where I45 and m are in units of 
10 45 g cm 2 and Mq , respective ly; the expected rang e of this 
quantity is between 0.5 and 2 JStella fc Vietrilll998ft . 

Recently, a model explaining the kilohertz QPOs 
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Figure 4. Left panel: The sideband separation frequency plot- 
ted vs. the lower kilohertz QPO frequency, Avsb- The full line 
is the quadratic relation between Avsb and the lower kilohertz 
QPO assuming the lower kilohertz QPO reflects the Keplerian 
frequency and assuming the sideband is formed due to a modu- 
lation caused by Lense-Thirring precession at that radius. Right 
panel: same as the left panel but for the upper kilohertz QPO. 



has been proposed that invokes the existence of a non- 
linear resonance between the vertical and radial epicyclic 
frequencies in an accre tion disc around a neutron star 
ijAbramowicz et afl l20031. Some versions of this model also 
include a resonance of thes e frequencies with the spin fre- 
quen cy of the neutron star jKluzniak et aljl2004t iLee et ail 
2004). The resonance model has been introduced as a way 
to explain the existence of prefer red values in the distri- 
bution of QPO frequency ratios jAbramowicz et all 120031 : 
but see iBelloni et al.N2004r ) as well as the commensurabil- 
ity of the neutron star spin and the kilohertz QPO fre- 
quency separation. So far it has not been explored whether 
a resonance mechanism could also explain the sidebands to 
the kilohertz QPOs. For instance, it may be possible that 
the Lense-Thirring precession at the radius at which the 
main resonance between the vertical and radial epicyclic fre- 
quencies occurs produces a modulation of the amplitude of 
the kilohertz QPOs. In principle, this modulation should 
produce equally strong sidebands, and should also produce 
sidebands on both QPOs. Alternatively, since for a rotating 
neutron star the azimuthal (Keplerian) frequency is larger 
than the vertical frequency, the sideband could be due to a 
resonance involving the azimuthal frequency at a radius at 
which a resonance between the vertical and radial epicyclic 
frequencies produce the upper and lower kilohertz QPOs. 
For instance, in this scenario a resonance between the az- 
imuthal and radial frequencies would produce a QPO at a 
frequency slightly higher than that of the upper kilohertz 
QPO. The advantage of this mechanism compared to an 
amplitude modulation is that it could in principle explain 
the fact that the upper sideband is stronger than the lower 
sideband. However, it remains unclear how a QPO at a fre- 



Sidebands is 4U 1636-53 5 



quency ~ 55 Hz higher than that of the lower kilohertz QPO 
would appear in this scenario. 

Previously, 4U 1636-53 was the only source out of the 
three sources where a sideband to the lower kilohertz QPO 
that showed no low-frequency QPO ll Jonker et akl koOO). 
Now, with the additional RXTE data, we could identify a 
low-frequency QPO at a frequency of ~43 Hz, similar to the 
one f ound in 4U 1728-34 and 4U 1608-53 (see I Jonker et alJ 
2000). The probable reason why this low-frequency QPO 
was not observed before is that the source did not sam- 
ple the states where this QPO is present often. Much of 
the new data were obtained in the state showing the low- 
frequency QPO. From TableQit can be seen that the lower 
kilohertz QPO frequency has to be below 850 Hz for the 
low-frequency QPO to be detected. Comparing the observed 
freque ncy distributions of the data set used bv lJonker et alJ 
( 2000) (Figure El shaded area top panel) with that used in 
this work (Figure|^|overall grey line top panel) it is clear that 
besides an increase in the amount of data, the lower kilohertz 
QPO frequencies in the range 700-850 Hz are much better 
sampled in the new dataset. A likely explanation for the fact 
that the distribution of the observed lower kilohertz QPO 
frequencies is different during the new RXTE observations 
from that observed in the data set used by Ijonker et alJ 
(2000), is that the X-ray flux of 4U 1636-53 has been de- 
creasing over the last few years (see Fig. E| bottom panel), 
allowing many of the additional new observations to be ob- 
tained at flux levels and at lower kilohertz QPO frequencies 
which were sparsely observed before in 4U 1636-53 with 
RXTE 
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